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Abstract. In this paper we study the problem of computing a Kolchin char- 
acteristic set of a radical differential ideal. The central part of the article 
is the presentation of algorithms solving this problem in two principal cases: 
for ordinary differential polynomials and in the partial differential case. Our 
computations are mainly performed with respect to orderly rankings. We also 
discuss the usefulness of regular and characteristic decompositions of radical 
differential ideals. In the partial differential case we give an algorithm for 
computing characteristic sets in the special case of radical differential ideals 
satisfying the property of consistency. For this class of ideals we show how to 
deal with arbitrary differential rankings. 



1. Introduction 

This paper is devoted to study radical differential ideals and their characteristic 
sets. The concept of a characteristic set introduced by Ritt and Kolchin is one 
of the most important notions in differential algebra. The problem of computing 
a characteristic set of a radical differential ideal represented by a finite set of its 
generators is not completely solved yet especially in the partial differential case. 
In the case of ideals in rings of polynomials in a finite number of variables this 
problem was studied and completely solved by Gallo and Mishra [HI [3 |H] ■ It was 
also investigated by Aubry, Lazard and Moreno Maza j2]- 

So, it is very natural to study this problem in rings of differential polynomials. 
The most important contributions of this article are: 

• an algorithm for computing characteristic sets of arbitrary radical differen- 
tial ideals w.r.t. orderly rankings in the ordinary case; 

• an algorithm for computing characteristic sets of radical differential ideals 
satisfying the property of consistency in the partial differential case. We do 
this w.r.t. both orderly fSections 14. II and 14.2(1 and arbitrary f Sect ion I4.3|l 
rankings. 

We also conjecture a method of solving this problem in non-ordinary cases for 
arbitrary radical differential ideals (see Section EJ. 

We use other techniques and methods than those used by Gallo and Mishra. 
However, their algorithm for computing a characteristic set of an algebraic ideal 
plays an important role in Algorithm^and Algorithm|2]of this paper. The methods 
developed by Sadik JJj help us to obtain several bounds for characteristic sets w.r.t. 
different rankings. 
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In ^Hl Sadik constructed an algorithm for computing a Kolchin characteristic 
set of a radical differential ideal w.r.t. an elimination ranking in the ordinary case. 
One can also get some bounds for orderly rankings from 18 , but these bounds are 
bigger then ours (see Remark 0J. Thus, the results of this paper (Theorem [S] and 
Theorem IB)) are new in comparison with Sadik's ones, also because Sadik did not 
give any algorithm for the partial differential case but we do. 

Ten years ago a technique for effective and factorization-free computations in 
the radical differential ideal theory was developed by Boulier, Lazard, Ollivier and 
Petitot (see pCI and In [51^17 Hubert continued to develop this problem and 
introduced the notions of characterizable ideal and characteristic decomposition of 
a radical differential ideal. This decomposition of the ideal helps us to solve many 
problems concerning the system of differential equations associated with the ideal 
and to test the radical membership. 

It should be emphasized that a characteristic decomposition of a radical differ- 
ential ideal does not give us full information about the ideal. In some important 
cases a representation of this ideal by characteristic components cannot replace a 
representation of the ideal by its generators as a radical differential ideal. For ex- 
ample, at this moment one cannot check the inclusion of a radical differential ideal 
to another radical differential ideal knowing only a characteristic decomposition of 
the first one (see j^l^^^]). This problem is closely related to the well-known 
Ritt problem. In this case it is necessary to know generators of the ideal and 
characteristic decomposition is partially useless. 

Nevertheless, in this paper we show that a characteristic decomposition of a 
radical differential ideal yields a lot of information about the ideal. Indeed, charac- 
teristic decomposition allows us not only to test membership to this ideal but also 
to compute its characteristic set in Kolchin's sense. Hence, the main contribution of 
this paper can be also considered as another application of a characteristic decom- 
position. Such a decomposition tells us a lot about a system of partial differential 
equations and we study what else one can do using a characteristic decomposition. 

We also study the case of partial differential polynomials and give an algorithm 
(Algorithm 121) for computing a Kolchin characteristic set. Theorem |H1 in Section 13 
provides a theoretical basis for this. Radical differential ideals, for which we propose 
the algorithm, satisfy the property of consistency (see Definition^)) w.r.t. an orderly 
ranking. We will see that these ideals are better than arbitrary radical differential 
ideals in the computational sense. 

In summary, although a characteristic decomposition cannot replace the ideal in 
computational sense, this decomposition allows us to compute such an important 
subset of a radical differential ideal as its characteristic set. 

2. Preliminaries 

2.1. Basic definitions. Differential algebra deals with differential rings and fields. 
These are commutative rings with 1 and a basic set of differentiations A = {6i, . . . , (5„} 
on the ring. The case of A = {6} is called ordinary. If R is an ordinary differential 
ring and y G i? we denote S'^y by j/'^'^^. The ring of differential polynomials was 
introduced to deal with algebraic differential equations. 

Recent tutorials on constructive differential ideal theory are presented in j^ , 19 . 
We also use the Grobner bases technique discussed in detail in ^ . The definition of 
the ring of differential polynomials in I variables over a differential field k is given in 
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|ll[li;^ini)] . This ring is denoted by fc{yi, ...,?/;}. We consider the case of char k — Q 
only. We denote polynomials by f,g,h,... and use the notation /, J, P, Q for ideals. 

We need the notion of reduction for algorithmic computations. First, we intro- 
duce a ranking on the set of differential variables of fc{yi, . . . ,yi}. Construct the 
multiplicative monoid O = ((jj^Mj^ • • • 6^'^ , ki ^ 0). The ranking is a total ordering 
on the set {9yi} for each 6* e O and 1 ^ i ^ I satisfying the following conditions: 

(1) 0u ^ u, 

(2) u^v=^9u^ 0v. 

In later discussions we suppose that a ranking is fixed. 

Let M be a differential variable in k{yi, . . . , ?/;}, that is, u — Oyj for a differential 
operator 9 = 61^62^ ■ ■ ■ S^" G 8 and 1 ^ j ^ I. A ranking is said to be orderly iff 
ordit > ordf implies u > v for all differential variables u and v. A ranking >ei is 
called elimination iff yi >e/ yj implies 9iyi >ei ^2% for all 9i,92 G O. 

The highest ranked derivative 9yj appearing in a differential polynomial / e 
k{yi, ... ,yi}\k is called the leader of /. We denote the leader by Uf. Represent / 
as a univariate polynomial in uy: 

/ = Ifu] + aiu"}-^ + . . . + a„. 

The polynomial If is called the initial of /. 
Apply any ^ e A to /: 

df 

6f = ^r-Suf + SIfvJi + Saiu^r^ + . . . + 
duf J -I J 1 

The leading variable of 5f is 5uf and the initial of 5f is called the separant of /. 
We denote it by Sf. Note that for all G 0, 7^ 1, each 9f has the initial equal to 
Sf- 

Define the ranking on differential polynomials. We say that / > g iff u/ > Ug or 
in the case of m/ = Ug we have deg„^ / > deg„^ g. Let F C k{yi, . . . , yi} be a set of 
differential polynomials. For the differential and radical differential ideal generated 
by F in k{yi, ...,?//}, we use the notation [F] and {F}, respectively. 

We say that a differential polynomial / is partially reduced w.r.t. g iff no proper 
derivative of Ug appears in /. A differential polynomial / is reduced w.r.t. g 
iff / is partially reduced w.r.t. g and deg„ / < deg„ g. Consider any subset 
A C k{yi, ...,!/;}. We say that A is autoreduced iff A n fc = and each element of 
A is reduced w.r.t. all the others. Every autoreduced set is finite (see ^2 Chapter 
I, Section 9]). For autoreduced sets we use capital letters A, B, C, .... 

We denote the product of the initials and the separants of the elements of A by 
I A and Sa, respectively. Denote I a ■ Sa by Ha. Let £1 be a finite set of differential 
polynomials. Denote by S°° the multiplicative set containing 1 and generated by S. 
Let / be an ideal in a commutative ring R. Let / : S°° = {a € R\3s € S°° : sa G I}. 
If / is a differential ideal then / : S°° is also a differential ideal (see |lllll6ininil9j V 

If we want to enumerate the elements of A we write the following: A = Ai, A2, . . . , A,^ 
Let A = Ai, . . . , Ar and M = Bi, . . . , Bg he autoreduced sets. Let the elements of 
A and B be arranged in order of increasing rank. We say that A has lower rank 
than B iff there exists k ^ r,s such that rank Ai = rank for I ^ i < A: and 
rankAfc < rankS^, or if r > s and rank Ai = ranki?i for 1 ^ i ^ s. We say that 
rank A = rankB iff r = s and rankA^ = rankS^ for 1 ^ i ^ r. 
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Consider two differential polynomials / and g in R ~ k{yi, . . . ,yi}. Let / be the 
differential ideal in R generated by g. Applying a finite number of differentiations 
and pseudo-divisions one can compute a differential partial remainder fi and a 
differential remainder f2 of / w.r.t. g such that there exist s G S'^ and h G 
satisfying sf = fi and hf = /2 mod / with /i and /2 partially reduced and 
reduced w.r.t. g, respectively (see ^ for definitions and an algorithm for computing 
remainders). 

Let A be an autoreduced set in fc{yi, . . . , y;}. Consider the polynomial ring 
k[xi, . . . ,Xn] with xi, . .. ,Xn belonging to QY for Y = j/i, . . . Let U,V C 
be the sets of "leaders" and "non- leaders" appearing in the autore- 
duced set A, respectively. We denote k[xi,. . . , Xn] by and the leader of Ai 
by UAi or Ui for each 1 ^ i ^ p. 

The notion of a characteristic set in Kolchin's sense in characteristic zero is 
crucial in our further discussions. This was first used by Ritt for prime differential 
ideals, but Kolchin introduced characteristic sets for arbitrary differential ideals. 

Definition 1. 11, page 82] An autoreduced set of the lowest rank in an ideal / is 
called a characteristic set of /. 

We call these sets Kolchin characteristic sets in order to omit confusions with 
different notions, e.g., in [Sl^] characteristic sets are used in Kolchin's sense and 
in some other senses. 

As it is mentioned in ^2 Lemma 8, page 82], in characteristic zero A is a 
characteristic set of a proper differential ideal / iff each element of I reduces to zero 
w.r.t. A. Consider the definition of a characterizable radical differential ideal. 

Definition 2. [S| Definition 2.6] A radical differential ideal / in k{yi, . . . ,yi} is 
said to be characterizable iff there exists a characteristic set A of / in Kolchin's 
sense such that / = [A] : H"^. 

The following definition makes a bridge between differential and commutative 
algebra. Let f be a derivative in k{yi, . . . ,yi}. Ay is the set of the elements of A 
and their derivatives that have a leader ranking strictly lower than v. 

Definition 3. ^2 in.8] A is coherent iff whenever A, i? G A are such that ua and 
ub have a common derivative: v = tp'^J-A — 4>ub, then SsipA — Sa4>B G (A„) : H'^ . 

We emphasize that a characteristic set of a differential ideal is a coherent au- 
toreduced set (see [III 1131131111 )■ 

2.2. Important assertions. Consider several important results concerning radical 
differential ideals in rings of differential polynomials. The technique described in 
Pllllj helps us to cover some properties of these ideals. 

Theorem 1. |11[ III. 8, Lemma 5] Let A be a coherent autoreduced set in k{yi, . . . , yi} 
Suppose that a differential polynomial g is partially reduced w.r.t. A. Then g G [A] : 
tffgG{A):H^. 

Note that Theorem ^ is also known as Rosenfeld's lemma. 

Tiieorem 2. Theorem 3.2] Let A be an autoreduced set of k[V] [U] . // 1 ^ (A) : 

S'^ then any minimal prime of (A) : 5^ admits the set of non-leaders of A, V , as 
a transcendence basis. More specially, any characteristic set of a minimal prime of 
(A) : has the same set of leaders as A. 
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Theorem 3. Theorem 4.5] Let A be a coherent autoreduced set of R = k{yi, ...,?/;} 
such that 1 ^ [A] : . There is a one-to-one correspondence between the minimal 
primes of (A) : in k[V][U] and the essential prime components of [A] : in 
R. Assume C; is a characteristic set of a minimal prime of (A) : . Then Ci is 
the characteristic set of a single essential prime component of [A] : ( and vice 
versa). 

Lemma 1. Let A = Ai, . . . , Ap be an autoreduced set in the ring k[xi^ . . . , Xm] = R 
and a characteristic set of (A) : L^ . Suppose that a polynomial f = flmX™ + . . . + 
Oq € R is reducible to zero w.r.t. A and the indeterminate Xt does not appear in Ai 
for each 1 ^ i ^ p. Then Oj is reducible to zero w.r.t. A for all j ^ m. 

Proof. Since / is reducible to zero w.r.t. A, there exists / G such that 

p 

1=1 

ti . m q 

Let gi — ^ ^j^t fo'^ each 1 ^ i ^ p. Thus, we have / • ^ Okx'l — ^ dkXf with 
j=i k=a k=0 

dk £ [Ai, . . . , Ap). Hence, L ■ Oi £ (A) for each 1 ^ i ^ m, that is, G (A) : 
Since A is a characteristic set of (A) : L^ , we have that aU are reducible to zero 
w.r.t. A. □ 



3. The ordinary case 

3.1. Bounds for the orders of characteristic sets. Denote |5| Algorithm 7.1] 
by X" Decomposition. An input of this algorithm is a finite set of differential poly- 
nomials F and its output is a set € = Ci, . . . , C„ of characteristic sets of char- 
acterizable ideals [C^] : forming a characteristic decomposition of the radical 
differential ideal {F} in k{yi, . . . ,yi}: 

= [Ci] :iI^n...n[C„] 

We will show how this decomposition helps to compute a characteristic set of {F}. 

The main idea of Algorithm ^ is to move our problem into commutative alge- 
bra. In order to do this we need to know bounds for the orders of elements of a 
characteristic set of a radical differential ideal. 

Let R = k{yi, . . . ,yi} with A = {6}. So, we are in the ordinary case. Dif- 
ferential dimension of a prime differential ideal P is the maximal number q such 
that P n fc{yij, . . . ,yi^} = {0}. If / is a differential polynomial then ord/ de- 
notes the maximal order of differential variables appearing effectively in /. Let 
A = Ai, . . . ,Ap be an autoreduced set. Define the order of A by the following 
equality: ord A = ordAi + . . . + ord Ap. If a set C is characteristic of the ideal P 
w.r.t. an orderly ranking then by definition the order of the ideal P equals ordC. 

Denote by P{s) the elements of P of the order less than of equal to s. The set 
P{s) is a prime ideal in the correspondent polynomial ring. As it is proved in [T^ 
the dimension of P{s) is a polynomial in s for s ^ h = ord P. More precisely, 
dim P(s) = q{s + 1) + h, where q is the differential dimension of the ideal P. 
Moreover, q = I — p and p is the number of elements of a characteristic set of the 
ideal P w.r.t. an orderly ranking. Thus, the number p does not depend on an 
orderly ranking. 
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Remark 1 . Almost all recent results in the dimensional theory for prime differential 
ideals (the theory of differential dimension polynomials) are presented in 

Lemma 2. |17l Proposition 17] Consider a prime differential ideal P of differential 
dimension q and of order h. For every subset {j/ij, . . . of {yi, . . . ,yi} the 

ideal P contains a differential polynomial in the indeterminates {yi-^ , ■ ■ ■ , } with 
order less than or equal to h. 

Lemma 3. jl7[ Lemma 23] Let P be a prime differential ideal, then P admits a 

dC' 

characteristic set C — Ci, . . . ,C„ such that — y — r does not lie in P, where Vk is 
a non-leader and li^k — ord(Ci, v^) for a pair (i, k) in {1, . . . ,p} x {1, . . . ,1 — p}. 

A characteristic set of a prime differential ideal is not unique, e.g., consider the 
ideal [x] G k{x, y} and the elimination ranking with x > y. Then the set 
is a characteristic set of the ideal [x] for any n ^ 0. Hence, if we do not propose 
any restriction, we will not be able to guarantee that a characteristic set of a prime 
differential ideal has restrictions for the orders of its characteristic sets. 

In order to omit this problem, in |17j Sadik used irreducible characteristic set 
and proved in |17l Lemma 19] that any prime differential ideal has an irreducible 
characteristic set. 

The following result will give us an opportunity to find a characteristic set of a 
prime differential ideal with good restrictions for the orders of its elements w.r.t. 
any differential ranking. 

Lemma 4. A prime differential ideal P in fc{j/i, ...,?/;} admits a characteristic set 
C = Ci, . . . , Cp with the following properties: 

(1) It is irreducible, that is, 

(a) Ci is an irreducible polynomial in k{yi, . . . , yi}, 

(b) each C^+i is irreducible in the ring Quot{k[V]/{Ci, . . . , Ci) : I°°)[U], 
where 

• V is the set of all variables appearing in the polynomials Ci, . . . , Ci, 

• J?° is the multiplicative system generated by the initials of the 
polynomials Ci, . . . ,Ci, 

• U is the set of that variables from C^+i that are not in V . 

(2) Let be a differential variable of the maximal order s appearing in C. 

(s) 

Let also y). does not appear in Ci, . . . , Ci-i but does appear in Ci. Then 
5m - i P. 

Proof. This is just a combination of the proofs of Lemma 13 and [171 Lemma 19]. 
Indeed, let us take an irreducible characteristic set C given by |17[ Lemma 19]. Let 
y^i^ be a differential variable of the maximal order s appearing in C. Let also y^f^ 
does not appear in Ci, . . . , Ci_i but does appear in d. 

If Si^t G P. Then replace Ci by Si^t in C and apply the process producing 
irreducibility from [13 Lemma 19] to the new C. The polynomials Q+i, . . . ,Cp 
will not be affected, since we do not change the ideal over which we factor. Thus, 

(s) 

we can proceed by induction on the degree of t/j in d and then by the maximal 
number s. □ 
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Theorem 4. Let P be a prime differential ideal of order h in k{yi, ...,?//} and > 
be a differential ranking. Then there exists a characteristic set C — Ci , . . . , C„ of 
the ideal P w.r.t the ranking > such that the order in yt of each Ci does not exceed 
h for all I ^t^ I. 

Proof. The proof of \V7i Theorem 24] will be valid for our purpose if we modify it in 
the following way. Denote the set {yk \ yk is not a leader of any Cj, 1 ^ j ^ n} by 
*Jt. If for some 6 G Q and t, 1 ^ t ^ I the variable dyt is the leader of some Cj then 
ord(Cg, yt) ^ h for all 1 ^ g ^ n immediately by Lemma 12 since C is autoreduced 
and ord9 ^ h. If not, we just take a polynomial / G k{yc fl P of order not 
greater than h. We have / is reduced w.r.t. C. Contradiction. 

Let yt g or and C be a characteristic set given by Lemma0] Now, the main idea 
is to reduce the polynomial /-,- e k{yc given by Lemma|21w.r.t. C where 0yc 
is the leader of Cj for some 9 G Q. Suppose that for some j, 1 ^ j ^ n, we have 
ord(C-, , yt) > h. 

Note that ord{fj,yCj) > ord(Cj, yc^)- Indeed, we may choose fj with If. ^ P. 
Hence, if OTd{fj,yc ) < ord(Cj,yc ) then the ideal P would contain an element 
that is not reducible w.r.t. the characteristic set C, because the ideal P is prime. 
Contradiction. Let ord(/j, ycj) ^ ord(Cj, )• 

Our improvement is to choose 

C = arg max ord(Cj , yt) 

and then take C; G C of the lowest possible rank instead of using induction as 
Sadik did. Let Ui = diyi for some 9i G Q and Ui be the leader of Ci. We have 
s = ord(Ci, yt) > h and r/ = ord(/i, yt) ^ ord(Ci, yi) = rc, where = fiiyi,'^) = 

Let us reduce each term (coefficients Oj, initial If. and its leader yi^^^) of fi 
first by Ci. We need to differentiate Ci q times and get the remainder / where 
^ q ^ r f — rc . Remember that fi depends only on yi,^, and their derivatives. 
Hence, applying further steps of reduction to the terms of / w.r.t. other Cj we 
need to differentiate them less than q times if Cj G C or not greater than q times 
if Cj ^ C. Indeed, the set C is autoreduced and the variables to reduce can come 
just from derivatives of variables from Ci. Moreover, Ci has the smallest rank in 
C. 

That is why after we reduce all leaders of C from / we get the polynomial 
depending effectively on y[^~^'^'' and s + q ^ s. Its initial w.r.t. yj''^''' is equal to 



^ci • ■ ■ ■ • ^hl • -^ci • • ■ • • • //. • 1 



dC, 
dy[ 



where ii , . . . , z„ , ji , . . . , j„ e Z^o and //^ is the remainder of //, w.r.t. C in the 
case oi r f > rc ■ Moreover, in the case of r/ = rc we are in PTTl Lemma 20] because 
of our choice of Ci and immediately get the inequality ord(/i, yt) ^ ord(Ci,yt). 
Continue with the "differential" case. 

Remember that P is a prime ideal. Hence, /^^ • . . . • • • . . . • S*^" ^ P, 
because Icj and Scj ^ P for all j, I ^ j ^ n. Moreover, P = [C] : and C is a 
characteristic set of [C] : H^. Thus, according to Lemma ^ the polynomial 
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is reducible to zero w.r.t. C that contradicts either to Lemma 0] or the discussions 
in the previous paragraph, because for a prime differential ideal the fact that an 
element is reducible to zero w.r.t. a characteristic set is nothing else the element 
belongs to the ideal. □ 

Remark 2. Note that Theorem 01 is actually a generalization of Sadik's result |17l 
Theorem 24] that was proved just for elimination rankings. 

So, now we are ready to prove a final bound for a characteristic set of a radical 
differential ideal and to obtain an algorithm computing this set. 

Theorem 5. Let I be a radical differential ideal in fc{yi, . . . , y;}. Let I = n"=i • 
H^, be a characteristic decomposition w.r.t. an orderly ranking with Ci = C} , . . . , Cf' . 
Let h be the maximal order of 'Ci for 1 ^ i ^ rt. Then the lowest differentially au- 
toreduced subset of an algebraic characteristic set of the ideal 

n 

I'=f]{0iCi,ord9iu^,^h):Hg^ 

i=l 

is a characteristic set of I w.r.t. the orderly ranking. 

Proof. Let / = fXj=i Pj be a minimal prime decomposition and / have a charac- 
teristic set C = Ci, . . . , Cp w.r.t. the orderly ranking. Suppose that for some z, 
1 ^ i ^ p, we have ord(Ci, yj-) > h, where 9yj. = Ui is the leader of Ci and 9 £ O. 
We may suppose this because the ranking is orderly. Note that i ^ 2, because the 
leader of the differential polynomial Ci appears among the leaders of characteristic 
sets of the ideals Pj . If not then multiplying all the lowest elements of character- 
istic set of the ideals Pj we get a non-reducible w.r.t. C element of /, that is a 
contradiction. 

Let li be the initial of with li G Pi, . . . ,Pt and li ^ Pt+i, ■ . ■ ,Pa- Concentrate 
our attention at Pt+i, • ■ • , Pa- Denote a characteristic set of Pj w.r.t. the orderly 
ranking by Mj. We have li is not reducible to zero w.r.t. Mj for t + 1 ^ j ^ a. 
According to Lemma ^ there exists an element with the leader u'^ in each Mj such 
that Ui = 9iju'^ for some 9ij € & and all j, t + I ^ j ^ a. 

Consider a particular Pr for some r, t + I ^ r ^ a. Let, for simplicity, u^ = 
9sys for all 1 ^ s ^ p. Let Ci depend on derivatives of y = {j/jj , . . . , } C 
{yi, . . . , yi-i} — Y, its leader Ui — 9iyi, and some non-leaders. Introduce a "shifted" 
function ord' as follows. For h' — max ord{Ci,yj), we put 

ord' 9yj = ord 9 + h' ~ ord{Ci,yj). 

Hence, the ord'-orderly ranking >ord' appears on the set QY. That is, if ord' u > 
ord' V then u >ord' v for all u,v QY. Choose any rankings on the sets of deriva- 
tives 9(y \ Y) and @{yi+i, . . . , yi}. Let B = . . . , Dp^ be a characteristic set of 
Pr w.r.t. the ranking 

e(y \ Y) >ei BY >ei y^ > el Q{y^+l, ■ • • , yi} . 
According to Theorem ^ we have ovd{Dq,yj) ^ h for all g, 1 ^ g ^ r and j. 

The polynomial Ci is reducible w.r.t. D. Then, there exists s, 1 ^ s ^ pr, such 
that one can apply a step of reduction to Ci w.r.t. Ds. Hence, ord'uD^ ^ ord' of 
some variable from d . Remembering the definition of ord' and the fact that C is 
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a characteristic set of / w.r.t. the orderly ranking we get OTd{Ds,yj) ^ ord(Ci, j/j) 
for aU j, 1 ^ j ^ i — 1. We have two possible cases. 

The first one is that ord{Ds, yj) = ord{Ci,yj) for some 1 ^ j ^ i — 1. Then, 
either ord(Ci, yi) ^ h and we have nothing to prove, or ord(Ci, yi) = h + m > h. 
Let us take care of this case. We have 

ord' {Ds,yj) = ord'(Ci, y-,) = ord' {Ci,yk) 

for all 1 ^ fc ^ z — 1. Let ud^ — Vi'^"'^ ■ Since ord'u^i^ ^ ovd' {Ds,yj), we obtain 
that ord' = ord'(Ci, y^^), that is, hg ~ ord(Ci,yi^). 

Represent d = Id (vi^^"^^"^ + oi (^i''^™') + . . . + a„;. So, according to 
the previous paragraph one can apply an algebraic pseudo-reduction to Id and aq 
w.r.t. Dg for all g, 1 ^ 5 ^ n^. We obtain a new reduced Ci multiplied by Id^. 
Remember that r ^ t + I and, thus, both Id and Id^ do not lie in P,. so that after 
this step we obtain new Ci with the initial does not belong to Pr, because the ideal 
Pr is prime. 

Since Ci is reducible to zero w.r.t. D, we get a new Dg that can reduce the 
polynomial Ci. If is a derivative of the leader of this Dg, we are done in this 

case, since the differential polynomial Dg belongs to k{yi, ■ ■ ■ ,yi} and ordDs ^ h. 
Taking into account that Iq, is not reducible to zero and continuing this finite 
process of reduction we get another Dg in P,. with the property oid{Dg,yj) < 
ord(Ci, yj) for all j, 1 ^ j ^ i — 1. That is the second case. In this one we can put 

a 

Ci = Dsa,ai 
a=t+l 

where Dg^^a comes from Pa by the just described procedure. 

Finally, multiplying Ci by Id we obtain an element of / that is not reducible 
w.r.t. C. Contradiction. Thus, we know an upper bound for the orders of C;, 
I ^ i ^ p. Moreover, by Theorem |31 we have max ordBj max ordCj. By 

Theorem^ and Lemma ^ we obtain that 

[Q] : i?^ n k[9,y,, orde, ^ h] ^ (^^ , ord 0^^, ^ h) : H^^ 

and ordui ^ h for each i, 1 ^ i ^ n. Thus, / n k[9iyi, ord6'i ^ /i] = /'. 

Here we used the fact that the set of leaders of any characteristic set of the ideal 
P coincides with the one of C, if P is a minimal prime of [Ci] : H^. This is true 
due to Theorem[21and TheoremOl If P is a minimal prime of / then P is a minimal 
prime of some characteristic component [Ci] : of / and vice versa. Hence, our 
problem is purely commutative algebraic now. 

In order to get a characteristic set of the ideal / it is sufficient to compute an 
algebraic characteristic set C of /' w.r.t. the induced ranking on the variables of 
this polynomial ring and then find in C the lowest differentially autoreduced subset 
C. This set differentially reduces the ideal /' to zero and some characteristic set of 
/ lies in /'. Thus, C is a characteristic set of /. □ 

Remark 3. The bound obtained in Theorem [3 most probably is not true for partial 
derivatives. At least, the method of "shifted" rankings does not work in non- 
ordinary cases. This will be shown in Example 



10 



MARINA KONDRATIEVA AND ALEXEY OVCHINNIKOV 



Remark 4. One can also get another bound using Lemma 4.3]. It is going to be 
the bound from ThcorcmElmultiplied by the number of differential indeterminates. 
For example, in fc{yi, . . . ,yi} it will be equal to I ■ h that is greater than h. So, we 
have got a better bound using our method. 

Remark 5. By means of the proof of Theorem |S1 one can also complete the proof of 
|18l Theorem 4.4] showing how to get an element of the ideal / having elements of 
Pi, 1 ^ i ^ a such that the new set becomes autoreduced. This is not clear from 
Sadik's paper. 

3.2. Algorithm. As an immediate consequence of Theorem|Slwe have the following 
algorithm. Let some orderly ranking be fixed. 

Algorithm 1. Ordinary Characteristic Set Computation 

Input: a finite set F of ordinary differential polynomials. 
Output: characteristic set of {F} in Kolchin's sense. 

• Let £ — x-Decom position (F) and £ = Ci, . . . , C„ with Ci = C/, . . . , Cf' . 

• Let h = max ordC^. 

l^i^n 

• Compute V = fX^M^^i ^^rdelu^, s^h):H^^. 

• C := an algebraic characteristic set of the ideal /'. 

• Return the differentially autoreduced subset of C with the lowest rank. 

Remark 6. The last steps of Algorithm ^ can be performed by means of com- 
putations discussed in and IH|- More precisely, one can compute each 
li = {9l Cf , ord 6f Ufjj ^ h) : H^. using the Rabinovich trick and the elimination 
technique. Then, the intersection of the ideals /' = HlLi computed. 
The solutions to these two problems are presented in [3] . Finally, an algorithm for 
computing an algebraic characteristic set of /' is given in [HlIZIIHl- 

Remark 7. Algorithm x-Decomposition used in Algorithm ^ can be replaced by 
Algorithm Rosenfeld_Groebner presented in [11 and 5 and implemented in Maple. 

4. Partial differential case 

4.1. Theoretical bases. We are not in the ordinary case now. Introduce a special 
class of radical differential ideals that are good in the computational sense. 

Definition 4. We say that a radical differential ideal / satisfies the property of 
consistency iff there exists a characteristic set C C / such that 1 ^ [C] : . 

It is clear that any characterizable radical differential ideal satisfies the property 
of consistency. Moreover, it follows from Theorem[21and Theorem|21that any proper 
regular differential ideal (an ideal of the form [A] : for a coherent autoreduced 
set A) satisfies this property. Consider an example of non-regular radical differential 
ideal satisfying the property of consistency. 

Example 1. Let A = x{x — 1), xy, xz in k{x,y,z} with x < y < z. We have 
1 ^ [A] : H'^ . Consider the minimal prime decomposition: 

{A}== [x]C^[x-l,y,z]. 

We see that A is a characteristic set of {A}. Then the radical differential ideal 
{A} satisfies the property of consistency. Nevertheless, since minimal primes of 
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{A} have different sets of leader, {A} is not a regular ideal due to Theorem |3 and 
Theorem 01 

So, we are ready to prove Theorem El Let 0yi be a differential variable in 
k{yi, . . . ,yi}. Then, by definition its order equals ord6'. 

Theorem 6. Let I be a radical differential ideal in k{yi, . . . ,yi} satisfying the 
property of consistency and a characteristic set <C d I with 1 ^ [C] : . 

(1) Let U be the set of leaders of C and U' be the set of leaders of any charac- 
teristic decomposition of L . Then U C U' . 

(2) Let L = nr=iP^] • characteristic decomposition w.r.t. an orderly 
ranking with <Ci — , . . . , Cf ' . Let h be the maximal order of differen- 
tial variables appeared in the elements of Ci for 1 ^ i ^ n. Then the 
lowest differentially autoreduced subset of a characteristic set of the ideal 
fXi^ii^lCj ,OTd0lu^j ^ h) : is a characteristic set of I w.r.t. the 
orderly ranking. 

Proof. We have [C] C / C [C] : i?^. Consider the minimal prime decomposition 
I = ^T=lP^■ We have 

n 771 

J=f|[Q]:i7^ = fl^- 

i=l i=l 

Some components of the characteristic decomposition may appear to be unnec- 
essary. Let I = HiLiPj] ■ -^c° be a minimal characteristic decomposition, that is, 
I 7^ PliLi i^^j l'^i] ■ fo'" I ^ i ^ k. Let U' be the union of leaders of Ci for 
1 ^ i ^ k. If P and Pj are prime ideals for each 1 ^ j ^ t and P D 0^=1 Pi then 
P D Pi for some 1 ^ j ^ t (see Proposition 1.11]). Thus, if P is a minimal 
prime of / then P is a minimal prime of [Ci] : for some 1 ^ i ^ k. 

We obtain that the set of leaders of any characteristic set of P is equal to 
those of Ci by Theorem (21 and Theorem O . Hence, the union of leaders of char- 
acteristic sets of minimal primes of / is equal to U'. Include [C] : into 
a characteristic decomposition of /. For this purpose represent C as an out- 
put of Coherent-Autoreduced algorithm (see Algorithm 5.1]). Thus, we have 
I =[C]: H^n [B2] : n . . . n [Mr] : iJ^. Denote the set of leaders of C by U. 

Let P be a minimal prime of [C] : H^. Then P is a minimal prime of {C} (see 
[3 page 644]). Thus, P is a minimal prime of /, because {C} C / C [C] : H^. 
Since P is a minimal prime of [Ci] : for some 1 ^ i ^ fc, then the set of 
leaders of Ci is equal to U and U C U' . So, we know an upper bound for the 
order of a characteristic set of /. Due to Theorem ^ and Lemma we have [Ci] : 
iJg" n k[0,y„ orde^ h] = {OlCl .ovdOlu^, h) : H^^ and ordui < h for each 
1 ^ i ^ n. Thus, we obtained the result because the end of the proof is the same 
as in Theorem 121 □ 

The main contribution of TheoremjSlis that our computations are moved into the 
ring of commutative polynomials in a finite number of variables. This is a crucial 
point in Algorithmic 

Remark 8. We see that the set of leaders [/ of C is not only a subset of U'. We have 
U is equal to the set of leaders of some characteristic component in any characteristic 
decomposition of /. Thus, U is "concentrated" in some characteristic component. 
We call it the "localization" property. 
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Remark 9. The fact that the set of leaders of C is equal to that of some characteristic 
component holds true for any differential ranking. This follows from the proof of 
Theorem |S1 

Look at the following example of a non-regular radical differential ideal illustrat- 
ing Theorem El 

Example 2. Consider the radical differential ideal / = {xy} in k{x,y} w.r.t. any 
differential ranking such that x < y. We have the following decomposition: 

I^[x]n[y]. 

The ideal / is not prime. Moreover, it is not neither characterizable nor regular. 
Nevertheless, the ideal / satisfies the property of consistency, because one can 
consider its characteristic set C = xy in Kolchin's sense with 1 ^ [xy] : x°°. 

As Theorem El tells us the set of leaders y of the characteristic set C must equal 
to the set of leaders of some component of a characteristic decomposition of the 
ideal /. This is the case for its second component [y]. In order to compute C we 
do not need to differentiate anything. We just calculate the reduced Grobner basis 
of {x) n (y), that is equal to xy, and output this set as a characteristic set of the 
ideal /, since xy is a differentially autoreduced set. 

Note that the condition that an ideal satisfies the property of consistency cannot 
be omitted in Theorem (HI To support this fact we give Example 

4.2. Algorithm. In conclusion, we obtain the following algorithm. Let some or- 
derly ranking be fixed. 

Algorithm 2. Characteristic Set Computation 

Input: a finite set F of differential polynomials such that {F} satisfies the 
property of consistency. 

Output: characteristic set of {F} in Kolchin's sense. 

• Let C = x-Decomposition(i^) and £ = Ci, . . . , C„ with d — , . . . , Cf' . 

• Let h ~ max max ord C\ . 

• Compute r = n-Li(^'C'^ord0^Up, <ih):H^^. 

• C := an algebraic characteristic set of the ideal /'. 

• Return the differentially autoreduced subset of C with the lowest rank. 

The last steps of Algorithm j^l can be performed in the same way as done in 
Algorithm m 

Remark 10. Algorithm x-Decomposition used in Algorithm (21 can be replaced by 
Algorithm Rosenfeld_Groebner presented in ^ and and implemented in Maple. 

Note that h from Algorithmj^is less than h from Algorithm^in general. Hence, 
the ideals satisfying the property of consistency are better in the computational 
sense than radical differential ideals in general. 

4.3. Arbitrary differential rankings. TheoremjSlcan be easily applied to com- 
putation of Kolchin characteristic sets of radical differential ideals satisfying the 
property of consistency not only in the orderly case. The statement of the theo- 
rem remains the same eliminating the word "orderly" . But in order to move the 
problem to a commutative ring in a finite number of variables we also need to use 
a characteristic decomposition w.r.t. an orderly ranking. 
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So, the algorithms coming from the proof of Theorem is the following. Let 
some differential ranking > be fixed. 

Algorithm 3. Characteristic Set Computation 

Input: a finite set F of differential polynomials such that {F} satisfies the 
property of consistency. 

Output: characteristic set of {F} in Kolchin's sense. 

• Let !8 = x-Decomposition(F) w.r.t. the ranking > and 05 = Mi,...,Mk 
with M, ^Bl,...,Bf\ 

• Let h = max max ord i3f . 

• Let £ — x-Decomposition(i^) w.r.t. an orderly ranking and £ = Ci, . . . , C„ 
with 'Ci = Cl,...,Cf . 

• Compute r = {^l^^{OlCi,oiAelu^, ^h):H^^. 

• C :— an algebraic characteristic set of I' w.r.t. the ranking induced by > . 

• Return the differentially autoreduced subset of C with the lowest rank. 

Corollary 1. There exists an algorithm for computing a characteristic set of a rad- 
ical differential ideal satisfying the property of consistency in the partial differential 
case. 

5. Examples 

We show how to apply Algorithmic] to a particular radical differential ideal. 

Example 3. Let I = {{x - t)x' , x'y', {x ~ t){z' + y')} in Q{t){x,y,z} with t' = I 
and an orderly ranking x < y < z. We have the following decomposition: 

I =[x~t,y']r\[x',z' + y']. 

So, the maximal order of variables appeared in this decomposition is equal to 
1. Hence, we need to compute the reduced Grobner basis G of the ideal /' = 
{x — t,x' — l,y') n {x' , z' + y'). This can be done by the elimination technique: 
G equals the intersection of the reduced Grobner basis w.r.t. the lexicographic 
ordering x < x' < y' < z' < w of the ideal {w{x — t), w{x' — I), wy', (1 — w)x\ 
(I - w){z' + y')) and the ring Q{t)[x , x' , y' , z']. 

Finally G = {x-t)x', x'{x' -1), x'y', {x ^t){z' + y'), x' z' - {z' + y'), y'{z' + y'). 
Then a characteristic set of /' equals C = {x — t)x', {x — t){z' + y') and by Theorem|H| 
a characteristic set of the radical differential ideal / is also C. 

The following example shows the difference between radical differential ideals 
with the consistency property and an arbitrary radical differential ideal. This can 
be considered as a "counter-example" for Remark [S] and Remark |21 

Example 4. Let A = x{x ~ 1), xy, {x — l)z. We have 1 G [A] : . Consider the 
minimal prime decomposition: 

{A} = [x, z]r\[x - 1, y]. 

We see A is a characteristic set of {A} with the set of leaders equals U = x, y, z. 
Thus, U does not necessarily correspond to the unique characteristic component. 
In this example U d x, z U x,y and the localization property is not valid. 

Note that Theorem is true for Example 01 in this case we have a restriction to 
the orders of the elements of a characteristic set A. Consider a "counter-example" 
for Theorem ini 
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Example 5. Consider a radical differential ideal defined by its characteristic decom- 
position: 

/=[x-l,y]n[:r,y("),z(™)+z/] 

in k{x, y, z} with x < y < z, an orderly ranking, and n ^ m. Both of these compo- 
nents are prime differential ideals, because they are generated by linear differential 
polynomials. In addition, since they are prime, these radical differential ideals are 
also characterizable (see [HI page 646]). One can show that a characteristic set of / 
is C = x{x — f), xy, (x — l)z(™+"). The radical differential ideal / does not satisfy 
the property of consistency and ThcorcmElis not true for /. Indeed, for m,n > 
we have m + n > max{m, n}. 

So, we see that the upper bound established in Theorem is wrong for some 
radical differential ideals not satisfying the property of consistency. Neverthe- 
less, Example 13 is in the ordinary case and Theorem can be applied. We have 
ord{x,?/("\ z^™) + y} ^ m + n and the maximal order of elements of C does not 
exceed m + n. 

Remark II. We see that the bound for radical differential ideals satisfying the 
property of consistency is lower than one for arbitrary radical differential ideals in 
general. That is why radical differential ideals satisfying the property of consistency 
are better in the computational sense. Indeed, Algorithm |21 works faster than 
Algorithm ^ since we are in a lower algebraic dimension in Algorithm [5] than in 
Algorithm m 

6. Conjecture 

The following example shows that the method of "shifted" rankings used in 
Theorem does not work in the partial derivative case. 

Example 6. Consider k{u} with A = {d/dx,d/dy,d/dz}. Let 

P \^^yyyy ^ ^-^xxyy ^ ^-^xxxx ~^ ^z: '^xxxxx ^-^xxy ^ ^xyy '^xxyyy] • 
Consider / — Uxxy + ^xyy + Uxxxxx + Uyyyyy + Uxxxxy + e P. One can prove 

that this polynomial cannot be an element of any characteristic set of P w.r.t. 
"shifted" rankings. However, / can be an element of a characteristic set of a 
radical differential ideal for which the ideal P is a minimal prime component. 

Remark 12. Example El is not a counter-example for Theorem El in non-ordinary 
cases. It is just a counter-example for the part of the proof of the theorem. 

Corollary 2. One need to develop another technique to get a bound for the order 
of a characteristic set of a radical differential ideal in partial derivative cases and 
radical differential ideals not satisfying the property of consistency. 

Nevertheless, we hope that this problem can be solved by means of increasing 
estimates of bounds for characteristic sets of radical differential ideals. So, we 
conjecture that if |A| = m then the order of each element of a characteristic set 
of a radical differential ideal / is bounded by q/i™, where h — max ordCi and 

l^^i^Cn 

/ = nr=i['^«] ■ ^ characteristic decomposition w.r.t. an orderly ranking and 

q is the number of different differential indeterminates t/j appearing as leaders in 
Ci for \ ^i ^n. 
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7. Conclusions 

We discussed an algorithm for computing a Kolchin characteristic set of an ar- 
bitrary ordinary radical differential ideal w.r.t. orderly rankings. In the partial 
differential case we also presented a solution to the problem of computing a charac- 
teristic set of a radical differential ideal satisfying the special property of consistency. 
These solutions are new and previously these problems were completely solved only 
in the non- differential case as to our knowledge. 

The authors hope that the technique obtained in this paper can be generalized 
to any radical differential ideals using the ideas we presented. Another natural 
way of generalizing these results is to complete our investigation of non-orderly 
rankings such as, for example, very important elimination ones in the case of partial 
differential polynomials. 
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